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ABSTRACT 
We report the discovery of a threshold in the Hi column density of Galactic gas clouds below which the 
formation of the cold phase of Hi is inhibited. This threshold is at Nhi — 2 y. 10^" per cm^; sightlines with 
lower Hi column densities have high spin temperatures (median T^ ^ 1800 K), indicating low fractions of 
the cold neutral medium (CNM), while sightlines with Nhi > 2 x 10^° per cm^ have low spin temperatures 
(median T^ ^ 240 K), implying high CNM fractions. The threshold for CNM formation is likely to arise due 
to inefficient self-shielding against ultraviolet photons at lower Hi column densities. The threshold is similar 
to the defining column density of a damped Lyman-a absorber; this indicates a physical difference between 
damped and sub-damped Lyman-a systems, with the latter class of absorbers containing predominantly warm 
gas. 
Subject headings: ISM; clouds — quasars; absorption lines — galaxies; high-redshift 
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1. INTRODUCTION 

The diffuse interstellar medium (ISM) contains gas over a 
wide range of densities, temperatures and ionization states. 
These can be broadly sub-divided into the molecular phase 
(e.g. Snow & McCall 2006) . the neutral atomic phase (e.g. 
iKulkarni & HeilesI 19881). and the warm and hot ionized 
phases (e.g. iHaffner et alj|2009l) . The neutral atomic medium 
(mostly neutral hydrogen. Hi) is further usually sub-divided 
into "cold" and "warm" phases (the "CNM" and "WNM", 
respectively). Typical CNM temperatures and densities are 
^ 40 — 200 K and > 10 cm^^^, respectively, with corre- 
sponding WNM values of > 5000 K and -- 0.1 - 1 cm"^. 
This was originally an observational definition, to distinguish 
between phases producing strong narrow absorption lines to- 
wards background radio-loud quasars and smooth broad emis- 
sion lines ( Clarklll965l) . Later, this separation into cold and 
warm phases was found to arise naturally in the context of 
models in which the atomic and ionized phases are in pres- 
sure equilibrium. Atomic gas at intermediate temperatures 
was found to be unstable to perturbations, and to migrate to 
either the cold or warm stable phases, given sufficient time 
to attain pressure equilibrium (Field et al. 1969; Wolfire et aT] 
11995). 

Radio studies of the Hl-21cm transition have played a vital 
role in understanding physical conditions in the neutral ISM. 
The Hl-21cm excitation temperature of a gas cloud (the "spin 
temperature", Tg) can be shown to depend on the local kinetic 
temperature (Tt-; lFieldlll959HLiszlil2001[). with T. x Tk in 
the CNM and T, < Tfc in the WNM (|Liszt||200lb. Since 
the Hl-21cm absorption opacity for a fixed Hi column density 
varies inversely with Ts, while the emissivity is independent 
of Tg (in the low opacity limit), Hl-21cm absorption stud- 
ies against background radio sources are primarily sensitive 
to the presence of CNM along the sightline, while Hl-21cm 
emission studies are sensitive to both warm and cold Hi. In- 
deed, the original evidence for two temperature phases in the 
neutral gas stemmed from a comparison between Hl-2 1 cm ab- 
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sorption and emission spectra (Clark 1965). 

Physical conditions in the neutral gas are also known to de- 
pend on the gas column density. At the low Hi column den- 
sities of the inter-galactic medium, Nhi < 10^^ cm^^, the 
gas is optically thin to ionizing ultraviolet (UV) photons and 
is hence mostly ionized, giving rise to the so-called Lyman-a 
forest (Rauch 1998). At higher Hi column densities, 10^^ < 
Nhi ^ 10^*^ cm~^, the Hi is optically-thick to ionizing pho- 
tons at the Lyman-limit, and the core of the Lyman-a absorp- 
tion line is saturated. Such "Lyman-limit systems" are par- 
tially ionized and typically arise for sightlines through the out- 
skirts of galaxies (Bergeron & Boisse 1991). At still higher 
Hi column densities, Nhi > 2 x 10^*^ cm^^, the Lyman-a line 
becomes optically-thick in its naturally-broadened wings and 
acquires a Lorentzian shape; such systems are called damped 
Lyman-a absorbers (DLAs; IWolfe et al]l2005 ). Damped ab- 
sorption is ubiquitous on sightlines through galaxy disks and, 
at high redshifts, has long been used as the signature of the 
presence of an in tervening galaxy along a quasar sightline 
( IWolfe et al.lll986l) . Finally, the molecular hydrogen fraction 
shows a steep transition from very low values t o greater than 
5% a t Nht ^ 5 X 10^° c m~^ in the Galaxy dSavage et alj 
Il977t iGillmon et all 12006)). Neutral gas is Ukely to be- 
come predominantl y molecu lar at higher Hi column densities, 
Nhi > 10^2 cm-2 (fSchavell2001l:lKrumholz et alJl2009h . 

While it is well known that a threshold column density 
(Nhi '^ 5 X 10^*^ cm^^) is required to form the molecular 
phase, for self-shielding against UV photon s in the H2 Ly- 
man band jSt echer & Williams 1967; Hollenb achet alJll971l; 
iFederman et alJll979D . the conditions for CNM formation are 
less clear. In this Letter, we report results from Hl-21cm ab- 
sorption studies of a large sample of compact radio sources 
that indicate the presence of a similar, albeit lower. Hi col- 
umn density threshold for CNM formation in the Milky Way. 

2. THE SAMPLE 

We have used the Westerbork Synthesis Radio Telescope 
(WSRT, 23 sources), the Giant Metrewave Radio Telescope 
(GMRT, 10 sources) and the Australia Telescope Compact 
Array (ATCA, 2 sources) to carry out sensitive, high spec- 
tral resolution (^ 0.26 — 0.52 km/s) Hl-21cm absorption 
spectroscopy towards 35 compact radio-loud quasars. Most 
of the target sources were selected to be B -array calibrators 
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for the Very Large Array, and have angular sizes < 5". De- 
tails of the sample selection, obs e rvations and data analysi s 
are given in iKanekar et all (l2003h . iBraun & Kanekari (l2005h . 
Roy et al. (2011, in prep.) and Kanekar & Braun (2011, 
in prep.). The final optical-depth spectra have root-mean- 
square (RMS) noise values of trms ~ 0.0002 - 0.0013 per 
1 km/s channel, with a median RMS noise of ^ 5 x 10^^ 
per 1 km/s channel. These RMS noise values are at off-line 
channels, away from Hl-21cm emission that increases the sys- 
tem temperature, and hence, the RMS noise. A careful band- 
pass calibration procedure was used to ensure a high spec- 
tral dyna mic range and excellent s ensitivity to wide absorp- 
tion [see IKanekar et all (12003.) and iBraun & Kanekari (120051) 
for details]. These are among the deepest Hi-21cm absorp- 
tion spectra e ver obtained (e.g. 'Dwa rakanath et al.l 120021: 
IKanekar et al.ll20 03; Begum et al. 2010) and constitute by far 
the largest sample of absorption spectra of this sensitivity. The 
use of interferometry also implies that the spectra are not con- 
taminated by Hl-21cm emission within the beam, unhke the 
case with si ngle-dish absorption spectra (IKanekar et al.ll20()3l; 
iHeiles & Tr oland 2003). 

Galactic Hl-21cm absorption was detected against ev- 
ery source but one, B0438-436. We used the Leiden- 
Argentine-Bonn (LAB) survey (iKalberla et a l. 20050 to es- 
timate the "apparent" Hi column density (uncorrected for 



TABLE 1 

The full sample of 35 sources. 



self-absorption, i.e. Nri = 1.823 x 10 
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dWilson et al.ll2009l) . where Tb is the brightness temperature) 
at a location adjacent to each target. The observational results 
are summarized in Table [Tj whose columns contain (1) the 
quasar name, (2) its Galactic latitude, (3) the Hi column den- 
sity (uncorrected for self-absorption), in units of 10^° cm~^, 
from the LAB survey, (4) the integrated Hl-21cm optical 
depth / rdV, in km/s, and (5) the harmonic-mean spin tem- 
perature Ts in K, estimated assuming the low optical depth 
Umit (T, = / TsdV/ / rdV). 

Finally, for the one source with undetected Hl-21cm ab- 
sorption (B0438-436), we quote a 3cr upper limit on the inte- 
grated Hl-21cm optical depth, assuming a Gaussian line pro- 
file with a full width at half maximum of 20 km/s, and the 
corresponding 3cr lower limit on the spin temperature. 

3. RESULTS: AN Nhi THRESHOLD FOR CNM FORMATION 

The left panel of Figure [T] shows the integrated Hl-21cm 
optical depth plotted against Hi column density, on a log- 
arithmic scale. It is clear from the figure that the inte- 
grated Hl-21cm optical depth has a rough power-law de- 
pendence on Nhi for Nhi > 2 x 10^" cm^^. However, 
there is a steep decline in / rdV for Hi column densities 
lo\ yer than this value. These p roperties were noted previously 
by IBraun & WalterbosI ( 119921) for individual spectral features 
(see their Figs. 3-6) rather than the line-of-sight integral. A 
similar pattern is visible in the right panel of the figure, which 
plots the spin temperature T^ against Nhi. 24 out of 25 
sightHnes with Nhi > 2 x 10^° cm^^ g^ow T^ < 600 K, 
while all ten sightHnes with Nhi < 2 x 10^° cm~^ have 
Ts > 550 K. The median spin temperature for sightlines 
with Nhi > 2 X 10^" cm^^ jg _ 240 K, while that for 
sightHnes with Nhi < 2 x 10^" K is - 1800 K. Thus, 
there appears to be a physical difference between sightlines 
with Hi columns lower and higher than the column density 
Nlim = 2 X 10^" cm^^. Note that the use of apparent Nhi 

^ http://www.astro.um-bonn.de/webaiub/english/tools_labsurvey.php 
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B0023-263 


-84.17 


1.641 ± 0.022 


0.025 ± 0.005 


3546 ± 694 


BOl 14-211 


-81.47 


1.380 ±0.034 


0.135 ±0.005 


559 ± 26 


B0117-156 


-76.42 


1.418 ± 0.034 


0.031 ± 0.004 


2538 ± 353 


BO 134+329 


-28.72 


4.272 ±0.029 


0.443 ± 0.002 


530 ±4 


B0202+149 


-44.04 


4.809 ± 0.027 


0.747 ±0.005 


353 ±3 


B0237-233 


-65.13 


2.078 ± 0.038 


0.294 ± 0.004 


388 ±9 


B0316+162 


-33.60 


9.431 ± 0.034 


2.964 ± 0.004 


175 ±1 


B03 16+413 


-13.26 


13.234 ±0.028 


1.941 ±0.003 


374 ±1 


B0355+508 


-1.60 


74.344 ± 0.079 


45.820 ±1.120 


89 ±2 


B0404+768 


18.33 


10.879 ± 0.032 


1.945 ±0.005 


307 ±1 


B0407-658 


-40.88 


3.363 ±0.018 


0.548 ± 0.007 


337 ±5 


B0429+415 


-4.34 


37.067 ±0.041 


10.879 ±0.007 


187 ±1 


B0438-436 


-41.56 


1.380 ±0.027 


< 0.020 


>3785 


B0518+165 


-11.34 


20.616 ± 0.035 


6.241 ± 0.007 


181 ±1 


B0531+194 


-7.11 


26.609 ± 0.036 


4.062 ±0.005 


359 ±1 


B0538+498 


10.30 


19.542 ±0.029 


5.618 ±0.003 


191 ±1 


B0831+557 


36.56 


4.469 ± 0.032 


0.483 ± 0.006 


507 ±7 


B0834-196 


12.57 


7.125 ±0.036 


0.973 ± 0.005 


402 ±3 


B0906+430 


42.84 


1.251 ±0.036 


0.051 ±0.003 


1342 ± 89 


B1151-348 


26.34 


7.732 ±0.031 


0.714 ±0.005 


594 ±5 


B1245-197 


42.88 


3.802 ± 0.037 


0.158 ±0.005 


1323 ± 40 


B1323+321 


81.05 


1.260 ±0.034 


0.083 ± 0.003 


836 ± 40 


B 1328+254 


80.99 


1.065 ±0.036 


0.021 ±0.003 


2729 ±411 


B 1328+307 


80.67 


1.197 ±0.030 


0.072 ± 0.002 


916 ±32 


B1345+125 


70.17 


1.957 ±0.024 


0.305 ± 0.005 


352 ±7 


B 161 1+343 


46.38 


1.318 ±0.034 


0.019 ±0.003 


3873 ± 565 


B 1641+399 


40.95 


1.044 ± 0.032 


0.009 ± 0.002 


6498 ± 1760 


B 18 14-637 


-20.76 


6.416 ±0.018 


0.997 ±0.007 


353 ±3 


B 1827-360 


-29.34 


8.163 ±0.018 


1.542 ±0.003 


290 ±1 


B 192 1-293 


-11.78 


7.312 ±0.017 


1.446 ±0.006 


277 ±1 


B2050+364 


-48.84 


27.638 ± 0.038 


3.024 ±0.010 


501 ±2 


B2200+420 


-5.12 


17.119 ±0.042 


3.567 ±0.016 


263 ±1 


B2203-188 


-51.16 


2.432 ±0.018 


0.248 ± 0.004 


537 ± 10 


B2223-052 


-10.44 


4.561 ±0.030 


1.034 ±0.003 


242 ±2 


B2348+643 


2.56 


70.522 ±0.051 


32.517 ±0.031 


119 ±1 



Notes: ' 



From the LAB emission survey, 
survey. 



From our Hl-21cm absorption 



only affects high-opacity sightlines, and has no significant ef- 
fect on our results. 

We used a number of non-parametric two-sample tests 
(e.g. the Gehan test, the Peto-Prentice test, the logrank test, 
etc), generalized for censored data, to determine whether 
the spin temperatures of sightlines with Nhi > Nlim and 
Nhi < Nlim are drawn from the same pare nt distribution. 
All te sts used the ASURV Rev. 1.2 package (Lavallev et alJ 
11992b which implements the methods of Feigelson & NelsoQ 
(t985). The use of multiple tests guards against any biases 
within a given test, resulting from the relatively small sam- 
ple size (10 system s with Nhi < Nlim) an d the presence 
of censored values (Feigelson & Nelsonlll985J) . Errors on in- 
dividual measurements were handled through a Monte-Carlo 
approach, using the measured values of T^ and Nhi (and the 
associated errors) for each sightline to generate 10'' sets of 
35 pairs of T^ and Nhi values. The statistical significance of 
each result (quoted below) is the average of values obtained 
in these 10^ trials for each test. The tests found clear ev- 
idence, with statistical significance between 4.2cr and 5.3cr 
that the spin temperatures for sightlines with Nhi > Nlim 
and Nhi < Nu m are drawn from d i fferen t parent popula- 
tions. Following iFeigelson & NelsonI (Il985h . our final result 
is based on the Peto-Prentice test, as this has been fo und to 
give the best results for very different sample sizes (Latta| 
Il981h . The two Tg sub-samples (with "low" and "high" Nhi) 
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Fig. 1 . — [A] Left panel: The integrated optical depth J rdV (in km/s), from our Hl-21cm absorption survey, plotted against apparent Hi column density 
(cm^^), from the LAB emission survey; the dashed cur'ves show the relation in eqn.|3]with {log{No), log{Noo), AV) = (20.0, 21.7, 20) and (20.3, 22.0, 10). 
[B] Right panel: The column-density-weighted harmonic mean spin temperature (K) plotted against apparent Hi column density (cm~^). 



are then found to be drawn from different parent distributions 
at 5.3cr significance; the probability of this arising by chance 
is - 2 X 10"^. 

It thus appears that sightlines with Nhi < Nlim — 1 x 
10^° cm^^ have systematically higher spin temperatures than 
sightlines with Nhi > Nlim- The spin temperature mea- 
sured here is the column-density-weighted harmonic mean of 
the spin temperatures of different "phases" along the line of 
sight (e.g. Kulkarni & Heiles 1988). For example, a sight- 
hne with Nhi equally divided between CNM and WNM, with 
T^ = 100 K and T^ = 8000 K respectively, would yield 
an average Ts ^ 200 K, while one with 90% of gas with 
T^ = 8000 K and 10% with T^ = 100 K would yield 
an average Ts — 900 K. In other words, a high T^ indi- 
cates the presence of a smaller fraction of the cold phase of 
Hl Thus, the fact that average spin temperatures are sig- 
nificantly higher on sightlines with low Hi column densities, 
Nhi < 2 X 10^° cm^^, indicates that such Iow-Nhi sight- 
lines contain low CNM fractions, far smaller than those on 
sightHnes with Nhi > 2 x 10^° cm~^. 

It is clear from Fig.fltA] that the integrated Hl-21cm opti- 
cal depth drops sharply at Nhi = Nlim, due to the decline in 
the CNM fraction at low Hi column densities. This can be ac- 
counted for by a simple physical model in which a minimum 
"shielding" Hi column density of WNM is needed for the for- 
mation of the cold phase in an Hi cloud. Assuming that the 
Hl-21cm optical depth in the WNM is negligible compared 
to that in the CNM, the equation of radiative transfer for a 
"sandwich" geometry can be written as 



Tb{^) = 



r,(l - e-^0 + ^e-^= (1) 



where the explicit velocity dependence of tw and tc has been 
omitted for clarity and the integrated Hi emission profile be- 



comes. 



ar. 



(2) 
(3) 



for a measured "apparent" column density, N^-^ (assuming 
negligible self-opacity), a threshold column density (where 
Tc — > 0), A^o ^ TiuTioAy, a saturation column density (where 
Tc — > oo), Noo ~ (Tc+Tujrm/2) Ay and an effective opacity, 
Tc- The effective opacity is related to the measured integrated 
opacity by the effective linewidth, AV^, as / rdV — t^AV^. 
The upper and lower dashed curves in Fig. [HA] show the 



above expression for (Nq 
10"' 

N. 



10 



20 



21 cm-2, AV = 20 km/s) and (Nq 



. N,. 
2 X 10^ 



5.0 X 



10 cm"-, AV — 10 km/s), respectively. Note that 
the effective linewidth only shifts the curves up and down in 
the figure. The fact that none of the data points of Fig. [11 A] 
lie to the left of the upper curve indicates that a minimum 
shielding Hi column of [Afo/2] ^ 5 x 10^^ cm^^ is needed 
for the formation of the cold neutral medium. Note that, in the 
above "sandwich" geometry, a total WNM column density of 
> lO'^'^ cm^-^ is needed before any CNM can be formed. 

4. DISCUSSION 

Most theoretical models of the diffuse interstellar medium 
treat it as a multi-phase medium, consisting of neutral and 
ionized phases in pressure equilibrium (e.g. F ield etaO 19691 : 
IMcKee & Ostriker 1977; Wolfire et al. 1995). In the McKee- 
Ostriker model, physical conditions in the ISM are regulated 
by supernova explosions, whose blast waves sweep up gas 
in the ISM into shells, leaving lar ge cavities full of hot ion- 
ized gas (McKe e & Ostrike^ll977l) . Cold atomic gas is then 
formed by the rapid cooling of the swept-up and shocked gas. 
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Fig. 2. — The column-density-weighted harmonic mean spin temperature 
(K) plotted against Galactic latitude. 

Soft X-rays from neighbouring hot ionized gas and stellar UV 
photons penetrate into the CNM and heat and partially ion- 
ize it, producing envelopes of warm neutral and warm ion- 
ized gas (WIM). Thus, a typical ISM "cloud" [see Fig. 1 of 
iMcKee & OstrikeJ (Il977l) 1 is expected to form at the edges 
of supernova remnants and super-bubbles containing the hot 
ionized medium (HIM). Such a cloud consists of a CNM core 
surrounded by a WNM envelope and a further WIM shell, all 
of which are in pressure equilibrium. The CNM core is almost 
entirely neutral, the WNM and WI M are partially ionized, an d 
the HIM is almost entirely ionized (iMcKee & Ostrikeril 19771) . 
A threshold Hi column density for CNM formation in an 
ISM cloud arises quite naturally due to the need for self- 
shielding against ionizing ultraviolet (and soft X-ray) photons 
(see also Schaye 2004 ). At low Hi columns, UV photons can 
penetrate all the way into a cloud core and heat (and partly 
ionize) the Hi. Only clouds that self-shield against the pene- 
tration of these high-energy photons can retain a stable CNM 
core. Self-shielding against UV photons is significant once 
the optical depth at the Lyman limit crosses unity, i.e. for 



N 



HI 



> 



and becomes more and more efficient 



with increasing Hi column density. Our results indicate that 
self-shielding only becomes entirely efficient at excluding UV 
photons from the interior of Hi clouds at a total Hi column 
density of Nhi ^ 2 x 10^° cm^^. Below this threshold, suf- 
ficient UV photons penetrate into the cloud interior to hin- 
der the survival of the cold phase. The simple "sandwich" 
model overlaid on Fig. [UA] suggests that a shielding column 
of ~ (0.5 — 0.75) X lO^'' cm^^ on each exposed surface of a 
cloud is sufficient to exclude UV photons, with the rest of the 
Hi then free to cool to lower temperatures. 

It is also possible that sightlines with low Nhi are sam- 
pling gas at higher distances from the Galactic plane, with 
lower metallicity and pressure. Two-phase models that in- 
corporate detailed balancing of heating and cooling rates ob- 
tain lower CNM fractions at lower pressures and metallic- 
ities ( fWolfire et a l.; 1995,). Observationally, .Kanekar et al.. 



( 120091) have found evidence for an anti-correlation between 
spin temperature and metallicity in high-z DLAs, with low- 
metallicity DLAs having higher spin temperatures, probably 
because th e paucity of metals yields fewer cooling routes 
(Kanekar & Chengalur 2001). If Iow-Nhi sightlines contain 
"clouds" with systematically lower metallicities (and/or pres- 
sures), these could have lower cooling rates and hence, lower 
CNM fractions. Unfortunately, we do not have estimates of 
either pressure or metallicity along our sightlines and hence 
cannot test this possibility. We note that all sightlines with 
high spin temperatures are at intermediate or high Galactic 
latitudes, 6 > 40° (see Fig. |2]l, suggesting larger distances 
from the plane. Conversely, the figure also shows that low 
spin temperatures (< 400 K) are obtained even for sightlines 
at high Galactic latitudes (6 ~ 40 — 70°). Measurements 
of the metallicity and pressure along the Iow-Nhi sightlines, 
through U V spectroscopy, would be of much interest. 

There have been earlier suggestions, based on semi- 
analytical or numerical treatments of self-shielding, that Hi 
clouds are pred ominantly neutr al for Nhi ^ lO^'' cm^^ (e.g. 
IViegas I 1995; W olfe et al.l l2005!). However, this is the first di- 
rect evidence for a change in physical conditions in Hi clouds 
at this column density. Note that this is a factor of several 
lower than the known threshold of Nhi ^ 5 x 10^" cm^^ 
for the formation of molecular hydrogen in Galactic clouds, 
set by the requirement of self-shielding against UV photon s 
at wavelengths in the H2 Lyman band ( ISavage et al.lll977h . 
Thus, there appear to be three column densities at which phase 
transitions occur in ISM clouds, at Nhi ^ 2 x 10^° cm~^ re- 
sulting in the formation of cold Hi, at Nhi ^ 5 x lO^'^ cm^^ 
resulting in the formation of molecular hydrogen, and finally, 
at Nhi > 10^^ cm^^, when most of the atomic gas is con- 
verted into the molecular phase. 

In this context, it is vital to appreciate that the abscissa of 
Fig. [T] refers to apparent Hi column density under the as- 
sumption of negligible self-opacity of the Hi profile. The 
saturation column density, Noo, is an observational upper 
limit to j TsdV and does not represent a physical limit 
on Nhi. In fact, the widespread occurren ce of the Hi 
self-absorption phenomenon within the Galaxy dGibson et alj 
120051) and the detailed modeling of high resolution extragalac- 
tic Hi spectra suggests th at self-opacity is a common occur- 
rence dBraun et alj|2009l) which can disguise neutral column 
densities that reach Nhi ^ 10^'^ cm^^. 

The original definition of a DLA as an absorber with Nhi > 
2 X lO'^*^ cm^-^ was an observational one, based on the re- 
quirement that the damping wings of the Lyman-a line be de- 
tect able in low-r esolution optical spectra of moderate sensitiv- 
ity dWolfe et al.iri986) . With today's 10-m class optical tele- 
scopes, it is easy to detect the damping wings at significantly 
lower Hi column densities, ^ 10^^ cm~^(e.g. iPeroux et alj 



I2003h : such systems are referred to as "sub-DLAs". There has 
been much debate in the literature on whether or not DLAs 
and sub-DLAs should be treated as a single class of absorber 
and on their relative importance in contributing to the cos 
mic budget of neutral hydrogen and metals (e.g . [Peroux et 
2003; Wolfe etal. 2005; Kulkai-ni etaD Em^. jWolfe et 
( I2005D argue that DLAs and sub-DLAs are physically differ- 
ent, claiming that most of the Hi in sub-DLAs is ionized and 
at high temperature, while that in DLAs is mostly neutral; this 
is based on numerical estimates of self-shielding in DLAs and 
sub-DLAs against the high-z UV background (Viegas 1995). 
Our detection of an Hi column density threshold for CNM 
formation that matches the defining DLA column density in- 
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dicates that DLAs and sub-DLAs are indeed physically differ- 
ent types of absorbers, with sub-DLAs likely to have signifi- 
cantly lower CNM fractions than DLAs, at a given metallicity. 
In summary, we report the discovery of a threshold Hi col- 
umn density, Nhi ^ 2 x 10^° cm^^, for cold gas formation 
in Hi clouds in the interstellar medium. Above this thresh- 
old, the majority of Galactic sightlines have low spin tem- 
peratures, Ts < 500 K, with a median value of ^ 240 K. 
Below this threshold, typical sightlines have far higher spin 
temperatures, > 600 K, with a median value of ^ 1800 K. 
The threshold for CNM formation appears to arise naturally 
due to the need for self-shielding against ultraviolet photons, 
which penetrate into the cloud at lower Hi columns and heat 



and ionize the Hi, inhibiting the formation of the cold neutral 
medium. 
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